---
title: "erick-vieira-script-tf-ad-ufpe-2018"
output: html_document
---

```{r cars}
library(tidyverse)
library(lmtest)
library(car)
library(ggplot2)
library(dbplyr)
library(arm)
library(XML)
library(faraway)
library(psych)
library(sandwich)
```

#Apontando o diret�rio
``````{r}
setwd("C:/Users/DELL/Desktop/An�lise de Dados")
getwd()
```

# Importando dados
``````{r}
bd <- read.table("erick-vieira-bd-tf-ad-ufpe-2018.txt", h=T,colClasses = c(rep("numeric",4),"numeric","numeric","numeric"))
bd                            
View(bd)
```

#Vari�veis do banco
``````{r}
names(bd)
```

# Verificando vari�veis
``````{r}
str(bd)
```

# An�lise explorat�ria dos dados
``````{r}
summary(bd)
head(bd)
var(bd)
cov(bd) 
cor(bd) 
```

# Teste de hip�tese para o Coeficiente de Correla��o
``````{r}
cor.test(bd$ajuda,bd$membrobrbm)
```

# Executando Diagrama de Dispers�o  
``````{r}
plot(bd$membrobrbm,bd$ajuda)            
abline(bd$membrobrbm,bd$ajuda) 
```
  
# Gr�ficos descritivos das vari�veis
``````{r}
hist(bd$ajuda, main = 'Histograma Ajuda Externa Brasileira')
```
``````{r}
hist(bd$membrobrbm, main = 'Histograma Membrobrbm')
```
``````{r}
hist(bd$percapita, main = 'Histograma Per Capita')
```
``````{r}
hist(bd$viaof, main = 'Histograma Viaof')
```


#Modelagem
``````{r}
# Modelo: regress�o linear m�ltipla = VD e VI + Controles (Semi log)
lm <- lm(data = bd, log(ajuda) ~ factor(membrobrbm) + viaof + log(percapita))
```

# Visualizando coeficientes do modelo
``````{r}
summary(lm)
```

# Construindo tabela de an�lise da vari�ncia
``````{r}
anova(lm)
```

# Construindo intervalos de confian�a para os coeficientes da regress�o
``````{r}
confint(lm)
```

# Verificando demais pressupostos do modelo
``````{r}
# Checando normalidade dos Res�duos
# M�todo 1) Construindo histograma de res�duos
hist(residuals(lm))
```

# M�todo 2) Checando suposi��o de normalidade dos erros por gr�fico de "Probabilidade Normal de Res�duos"
``````{r}
windows()
qqnorm(residuals(lm), ylab="Res�duos",xlab="Quantis te�ricos",main="")
qqline(residuals(lm))
```

# Checando heterocedasticidade
``````{r}
plot(lm, which = 1) 
```

# Checando multicolinearidade
``````{r}
require(faraway)
vif(lm) 
```

# Corrigindo heterocedasticidade com estimador White-Huber
``````{r}
coeftest(lm,vcov. = vcovHC)
```


``````{r}
